function realu = U2(x,t)
        s = pi/2*ones(size(x));
        for k = 1:2:201
          s = s + 4*exp(-k^2*t)*cos(k*x)/k/k/pi;
        end
        realu = s;
end
